Mon. Not. R. Astron. Soc. 000, 000-000 (1997) 



Printed February 1, 2008 



(MN 1^1^ style file vl.4) 



Gravitational lensing and the angular-diameter distance 
relation 



Fedja Hadrovic and James Binney 

Theoretical Physics, 1 Keble Road, Oxford OXl 3NP 



' l^ebruary 1, 2008 

m 



> 
O 



ABSTRACT 

We show that the usual relation between redshift and angular-diameter distance can be 
derived by considering light from a source to be gravitationally lensed by material that 
lies in the telescope beam as it passes from source to observer through an otherwise 
empty universe. This derivation yields an equation for the dependence of angular- 
diameter on redshift in an inhomogeneous universe. We use this equation to model the 
distribution of angular-diameter distance for redshift z = 3 in a realistically clustered 
cosmology. This distribution is such that attempts to determine qo from angular- 
diameter distances will systematically underestimate go by 0.15, and large samples 
would be required to beat down the intrinsic dispersion in measured values of go- 
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INTRODUCTION 
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'■"iS'lie large-scale structure of the Universe is believed to be 
iclosely approximated by one of the Friedmann-Lemaitre cos- 
Oiological models. These are characterized by the values of 
^Tfl iree parameters: the matter density, the radius of curva- 
ture of spatial sections and the cosmological constant. De- 
J^erminatio: i of these values has long been considered one of 
^he fundar ST^aTTasS'Tr'oBservanonar^osmoTog^ 
JlJ^onnection a potentially key observable is the relationship 
r^etween redshift z and the angular-diameter distance D{z), 
j^hich is defined to be the ratio of the linear diameter of 
■ ■ an object to the angular diameter that it subtends when 
observed at reds hift z. There have been may attempts to 
determine D{z) (Bandage 1988; Kellermann 1993; Crawford 



Schneider, Ehlers & Falco, 199z), we obtain a much sim- 
pler analysis and an equation (^ that involves the cos- 
mic density field rather than the cosmic metric or potential. 
This simplicity enables us to evaluate the effects of inho- 
mogeneity for the case of realistic cl ustering, rath er than 
the case of either weak perturbations (Sasaki 1987) or ran 



domly distributed point m asses ( Zel'dovich, 1964f Refsdal 
197C|; |Dyer fc Roeder, 1974 [Watanabe fc Tomita 199q ; [Kan" 



towski, Vaughan & Branch 1995)' 



1995) which has recently prompted a wider discussion about 



the feasibility of the method (Nilsson et al. 1993 



Dabroski 



Lasenby fc Saunders 199f: ; Kantowski, Vaughan & Branch 
1995| ; pteph anas'ITSaha 1995[ ). The form of D(z) depends on 
the geometry of the Universe in the sense that the larger the 
curvature K is, the smaller D is at a given redshift. That is, 
the more positively curved the Universe is, the more slowly 
the angular size of an object decreases as it is moved away 
from the observer. 

In an inhomogeneous universe deviations of the met- 
ric from the Friedmann-Lemaitre form give rise to fluctu- 
ations in measured values of D at fixed z. Previous anal- 



yses of this effect (Zel'dovich 


, 1964, Refsdal, 1970l Dyer & 


Roeder, 1974 


Sasaki 1987 




r 


Watanabe & Tomita 199(]|) 



worked di rectly from general relativity and produced re 
suits of CO isiderable complexity. By delegating relativistic 
considerations to the theory of gravitational lensing ( e.g. 



Over the last decade there has been a growing awareness 
of the importance of gravitational lensing for observations 
of high-redshift objects. Gravitational lensing and the de- 
pendence of D on z are two sides of the same coin: both 
phenomena are caused by the tendency of matter that lies 
between the observer and a distant object to focus radia- 
tion from that object, thereby increasing its apparent size 
and brightness. In Section 2 we demonstrate this connec- 
tion quantitatively by showing that the standard formula for 
D{z) in a Friedmann-Lemaitre universe can be obtained by 
applying conventional lensing theory to a Universe in which 
there is matter only in the telescope beam towards the ob- 
ject under study. In Section 3 we describe our model of the 
clustered cosmic density field, and in Section 4 we use this 
model to calculate probability distributions for D{z) from 
objects of various linear sizes. Section 5 sums up. 

We throughout use the convention that D and D, re- 
spectively, denote angular-diameter distance before and af- 
ter lensing is taken into account. 

Since luminosity distance Dl is rigorously related to 
angular-d iameter distance by Dl = D/{1 + z)^ (Ethering- 



ton 1993), our distributions of values of D imply identical 



distributions of values of 
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Figure 1. A sequence of gravitational lenses 



The unit system we use is based on G — c — Ho ~ 1, 
which significantly simphfies the equations in cosmology. All 
lengths quoted are scaled to Hq — 100/ikms~^ Mpc~^. 



2 THE ANGULAR-DIAMETER DISTANCE 
RELATION FROM LENSING 

2.1 A sequence of gravitational lenses 

Fig. 1 shows a ray that is deflected through angles ai and 02 
by two lenses, which it passes at impact parameters ^1 and 
^2, respectively. From the figure it is immediately apparent 
that 

6 = S^|i-a(6)5i3-a(6)523. (1) 

The generalization of this equation to an arbitrary number 
of lenses is easily seen to be 



D ~ 



i = l 



(2) 



Light that passes at radius vector ^ through a disc 
of matter that has uniform surface density E is deflected 
through an angle 



(3) 



Hence in this case the impact parameters are all parallel and 
satisfy 



Di 



j-i 



(4) 



Finally, if is a diameter of an object, then the angular- 
diameter distance of that object is Di = {(,1/ D\), so on 
taking the modulus of and dividing through by (^i/-Di) 
we conclude that true angular-diameter distances Di satisfy 

Dj = Dj - 47r ^ A A, (5) 



2.2 Application to a homogeneous universe 

We now we show that equation reproduces the famil- 
iar angular-diameter distance equation for a Friedmann- 
Lemaitre universe. We consider an empty universe. In such 



a universe there is nothing to single out a unique rest frame 
at any given event, so redshift is not uniquely related to dis- 
tance. This permits us simply to adopt the relation s{z) be- 
tween proper distance and redshift in a Friedmann-Lemaitre 
universe. We have that a{z) satisfies [Scheneider et al., eq. 
(4.47b)] 



-1/2 



(6) 



ds 
dl 

From the gravitational-focusing equation [Schneider et 
al. eq. (3.64)] in the case of empty space (vanishing Ricci 
tensor and shear) we have 



d^D 

dr2 



= 0, 



(7) 



where r is an affine parameter for the light beam. In terms 
of the wavenumber, k, we have ds/dr cx fc oc 1 + 2:, so we 
may use equation (g) to eliminate r from (Q) in favour of s. 
We then find that the focusing equation states that in our 
empty universe, as a function of z, angular-diameter distance 
D satisfies 



d'D{y,z) 
dz'^ 



= 0. 



(8) 



We also have the initial condition [Schneider et al. eq. (4.53)] 

db{y,z 



dz 



{1 + z)-\i + ^zy 



-1/2 



(9) 



Now we fill the telescope beam with the normal matter 
density of a Friedmann-Lemaitre universe and use equation 
(P) to calculate the angular-diameter distance of an object 
at 'redshift' z. We first take the limit of equation (P) in which 
there are an infinite number of discs. Since in our units the 
current critical density is 3/(87r), the disc that lies between 
z + Az and z has surface density 



S= (I + 2) -^-T-Az, 

Stt az 



(10) 



where Q, is the usual density parameter. With equation ^ 
this becomes 



E(^) = 



30 1 + z 



■Az. 



(11) 



Stt ^1 + Q.z 

Inserting this expression for E into equation (^) and pro 
ceeding to the limit dz 



D{z)=D{z) 



we find 
1+y 



dy , 

"VT+TT^ 



D(y)D{y,z) 



(12) 



We now convert this integral equation for D{z) into a dif- 
ferential equation. Differentiating we find 



dD _ dD 

dz dz 

d^D _ d^D 

dz^ dz^ 

-1^ 



„ ''7f 

^ VI + iiz 



vi 



i^){piy,z) 



(13) 



Diy){—D{y,z) 



Combining these equations and taking advantage equa- 
tions (tel) and (M), we recover the standard equation for 
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Figure 2. Averaging of the matter distribution over the beam 
cross section 



the angular-diameter distance in a conventional Friedmann- 
Lemaitre universe: 

D BD 

{l + z){l + Qz)-^+{^nz+lQ + 3) — + ^QD = 0.(14) 



2.3 Application to an inhomogeneous universe 

The most important feature of the above derivation is that it 
does not depend on Q being constant. In the first two terms 
of equation ( p^ ) Q appears as a result of the reparametri- 
sation (s t-^ z). It is not related to the local matter dis- 
tribution and can be thought of as the averaged density 
parameter (Q). The parameter Q in the last term, fr^D 
is related to the local matter density and comes directly 
from the gravitational lensing calculation. Therefore, in the 
case of a locally inhomogeneous universe that approaches a 
Friedmann-Lemaitre model in the large-scale limit, we can 
write 

(1 + ^)(l + (Q) + (i {^) z+ i {n) +3)^ 



dz 

+ lQ.{z)D = 0. 



(15) 



Note that Q,{z) describes the comoving matter density be- 
cause the physical density is p{z) — ^f^(l + z)^. 

Simply replacing Q by (Q) in all but the last term of 
equation (|l4) does not allow for a complete discussion of 
the effects of inhomogeneity on images: in addition to be- 
ing magnified by matter within the beam, images will be 
distorted and may be even split into multiple images. We 
have neglected these potentially important effects by (i) as- 
suming that the material that lies between redshifts z + dz 
and 2 forms a uniform disc, and (ii) neglecting shear that 
is induced by clumps of material that lie outside the beam. 
Futumase & Sasaki (1989) and Watanabe & Sasaki (1990) 
show that as long as the scale of inhomogeneities is greater 
than, or equal to galactic scale, shear does not contribute 
significantly to focusing. 

By contrast, the assumption that the beam is filled by 
a series of uniform- density discs constitutes a non-trivial 
approximation about the matter distribution in the beam, 
namely that we may average the density across the beam as 
shown in Fig. |^. 



3 STATISTICAL MODEL OF THE FIELD n{r) 

We now investigate the predictions of the generalized 
diameter-distance equation (^|). For this investigation we 
require a statistical description of the density field along the 
telescope beam. This is a random field, which we think of 
as a function of comoving distance x. We assume that fl{x) 
follows a log-normal distribution - see Coles & Jones (1991) 
for a discussion of the characteristics and advantages of the 



log-normal distribution in cosmology. We confine ourselves 
to the case of a critical-density universe: (Q) = 1. With these 
assumptions Q{x) is given by 



(16) 



where e{x) is a Gaussian random field. Without loss of gen- 
erality we set {e} = 0. 

We define the two-point correlation function, 5/ of a 
field f{x) by 

{f{x' + x)f(x')}-{f{x')f 



C/(^) 



(17) 



The correlation functions of the fields fl{x) and e(a::) are 
related by 

exp{a^^^{x)) - 1 



^n{x) 



(18) 



where — (^e^'^ is the variance of the Gaussian field. 

The Gaussian field e{x) is determined by its power spec- 
trum Pe{k), which is essentially the Fourier transform of 5^: 

'' ■ (19) 



n(fc) = §^ J dxe'^^U^). 

Hence, if we know ^n{x), we may construct realizations of 
O, by determining £,e{x) from equation (|l^) and then using 
equation ( p^ ) to determine Ps{k). 

The galaxy correlati on function may be approximated 
by (Padmanabhan, 1993) 



C(0 = - 



(20) 

where 7 ~ 1.8 and the correlation length is rc — 5.5/i~^Mpc. 
The correlation function of the density field is often assumed 
to have the same form, but a different amplitude. The bias 
factor is introduced by setting 

^GALAXIES 



b = 



(21) 



^MATTER 

Measurements indicate that 1 ^ b ^ 2. Hence we require 
such that £,n{0) — 1 and 



We have adopted the form 



Cn(0 



(22) 



(23) 



with anro = 6 ^^^rc- For b = 1.5 we find anro = 
4.5/i"^Mpc. The meaning of ro will be discussed later, but 
we immediately see that for small ro the model approximates 
the divergent galaxy correlation function better. 
From equation (|l^) we find 



at \ rg + r2 
From (|l^ the power spectrum is 



(24) 



1 

2^ 
1 

k 



da; In — = — e 



exp(— fcro) — exp(— fcroe'^^''^) 



(25) 
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Angular diameter distance 
distribution at z=3.00 



.1 



0.5 1 1.5 

Figure 3. The distribution of angular-diameter distance at 2 = 3 

The significance of ro now emerges: it determines how 
quickly |P£(fc)| approaches zero at large k. The smaller the 
value of ro the larger must be the wavenumber fcmax up 
to which we must sum the discrete Fourier transform from 
which we obtain realizations of Q{x). Physically, we should 
think of ro as the scale on which the matter distribution 
is smoothed by the finite width of our telescope beam and 
the diameters of the objects we are looking at. If we take 
ro = 10h~^ kpc, we have = 12.21. 

Due to computational constraints and limitations on 
sampling imposed by Nyquist's theorem, it was impractica- 
ble to generate a single random field on the range < 2: < 3. 
Instead, we divided this interval into 100 subintervals and 
create a scaled random field on each of them. This procedure 
destroys correlations between different intervals but these 
are physically unimportant because the correlation function 
is negligible at such large distances. 



4 RESULTS AND DISCUSSION 

4.1 Distribution of angular-diameter distances 

Once a realization of n{x) has been constructed, it is 
straightforward to solve equation ( p^ ) for D at any given va- 
lue of z. We repeated this operation for approximately 4000 
realizations of Q{x) to determine the distribution of angular- 
diameter distances at 2 = 3. Fig. ^ shows this distribu- 
tion. The distances are rescaled to the standard Friedmann- 
Lemaitre value 



1 - 



1 



(26) 



An important point on the graph is the Dyer-Roeder dis- 
tance corresponding to an empty light beam: 



(l + z)5/2 



(27) 



which for 2 = 3 gives D°^/D^^ = 1.55. 

We see that the distribution is strongly peaked on the 
Dyer-Roeder side of D^^ , with a long tail on the Friedmann- 
Lemaitre side. This is expected because regions within which 



ro=10^h'Mpc 




ro=10^h"'Mpc 





ro = l h Mpc 



-..III 



I 



D/D 



D/D 



Figure 4. Effeet of smoothing on the angular diameter distance 
distribution 



the density is below average occupy the great majority of the 
volume of the Universe. Hence, many light paths sample only 
low-density regions and the distribution in Fig. ^ is shifted 
towards D^^ . However, when the light beam does encounter 
a galaxy or other matter aggregation, it is strongly lensed. 
These events decrease the diameter distance and give rise to 
the tail on the Friedmann-Lemaitre side. 

It is important to understand the impact that smooth- 
ing of the matter distribution has on our results. It is com- 
putationally convenient to investigate this for an unreal- 
istic case: we take the correlation length to be 100 times 
its true value. That is, we investigate the case in which 
""nro = 450/i~^Mpc. Fig. ^ shows our results. 

For large ro the matter distribution is rather homoge- 
neous, so the distribution of D is narrow and peaked near 
D^^. As ro is decreased the universe becomes strongly in- 
homogeneous and the distribution of D becomes broader. 
Simultaneously, its peak shifts towards the empty-beam dis- 
tance 



4.2 Implications for go measurements 

One of the most important undetermined quantities of cos- 
mology is go, the deceleration parameter. For a fiat universe 
(K = 0), the angular-diameter distance D is related to go 
by 

D{qo,z) ^ 



1+2 



where 



R 



du 



li (nu'J' + (1 + go - 3n/2)ii2 + in - go 



a/2' 



(28) 



(29) 
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Suppose we attempt to use ( psj ) to determine go from an 
observationally determined value of D{z). We Eissume that 
the true values of the cosmic constants are those with which 
we have been working: S7=1,A = A" = 0, and thus that the 
true value of go is go = 5- Putting qo ~ \ + Sqo in equation 
we find 

Substituting this value of R into (ji^) gives 

-^[^^ = 1-0.15 5go. (31) 

This equation relates the error, Sqo, in the inferred value of 
go to the ratio of the measured value of D to the value D^^ 
that it would have if the Universe were homogeneous. The 
distribution of D/D^^ shown in Fig. H is centred on 1.025 
and has spread ~ ±0.06. By equation ( pi| ) the error in go to 
which this gives rise is 

5go = -0.17 ±0.4. (32) 
In connection with this result three points should be made: 

• We see that the conventional method of determining go 
from the angular-diameter redshift relation provides a biased 
estimator of go that will return significant underestimates of 
the true value. 

• Even perfect measurements of D{z) will return values 
of go that are widely scattered. The breadth of this scatter 
is such that an accurate determination of go would require 
an extremely large sample and a sophisticated statistical 
analysis of the data. 

• The errors in go to which inhomogeneities give rise de- 
pend on the scale of observed objects because this scale 
determines the effective spectrum of the inhomogeneities. 
Larger objects will yield smaller errors. 

This last point is unfortunate because, as Kellermann 
(1993) has emphasized, small objects are much more likely 
to constitute standard measuring rods than large objects, 
such as giant radio sources, whose linear sizes are likely to 
be sensitive to the mean cosmic density. 



5 CONCLUSION 

We have used the theory of gravitational lensing to derive 
the conventional relation between angular-diameter distance 
and redshift in a Friedmann-Lemaltre universe. The value 
of this derivation is that it is simple and shows that the ten- 
dency of the angular diameter of a distant object to increase 
with Q. arises because rays coming from the object are fo- 
cused by matter that lies within the telescope beam. Hence, 
the angular diameter of an object is sensitive to the precise 
disposition of matter in the neighborhood of the telescope 
beam: move matter just out of the beam and the apparent 
size of the object will diminish. Equation (^) expresses this 
fact mathematically. 

Since the Universe is strongly inhomogeneous on small 
scales, telescope beams to different objects at the same red- 
shift will contain significantly varying quantities of matter, 
and the apparent diameters of physically identical objects at 



a common redshift will vary. This variation gives rise to scat- 
ter in the angular-diameter distances D of a set of objects 
that lie at a common redshift. 

We have modelled the distribution of the values of D 
of objects at redshift z = 3 by assuming that the cosmic 
density field follows a lognormal distribution that matches 
the observed clustering of galaxies for bias parameter b = 
1.5. The distribution of D is very skew, with its peak at a 
value that exceeds that associated with the corresponding 
homogeneous universe, D^^, and a long tail to values smaller 
thauD^^. In consequence of this skewness, the conventional 
technique for measuring go from measurements of D{z) will 
systematically underestimate go. 

The width of the distribution of D at given z depends 
upon the assumed power spectrum P{k) of the cosmic den- 
sity field. The true power spectrum is thought to have con- 
siderable power on small scales, and this power will generate 
a very broad distribution of D for objects of small angular 
size. When the angular diameters of highly extended ob- 
jects are measured, only power on scales comparable to or 
larger than the linear size ro of the objects will contribute 
to the scatter in D. Ifence such measurements will yield less 
scattered values of D. For ro — 10/i~^kpc we estimate that 
D will scatter by ~ ±6% at redshift 2 = 3. Unfortunately, 
even this small scatter will cause the derived values of go to 
scatter by as much as ±0.4. The scatter in values of go that 
are derived from angular-diameter distances to parsec-sized 
objects such as those studied by Kellermann (1993), will be 
very much larger still. 
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